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ABSTRACT 

We describe modifications to the joint stepwise maximum likelihood method of Cole (2011) 
in order to simultaneously fit the GAMA-11 galaxy luminosity function (LF), corrected for ra¬ 
dial density variations, and its evolution with redshift. The whole sample is reasonably well-fit 
with luminosity {Qe) and density (Pg) evolution parameters Qe,Pe ~ 1-0,1-0 but with signif¬ 
icant degeneracies characterized by Qe « 1.4 — O.dPg. Blue galaxies exhibit larger luminosity 
density evolution than red galaxies, as expected. We present the evolution-corrected r-band 
LF for the whole sample and for blue and red sub-samples, using both Petrosian and Sersic 
magnitudes. Petrosian magnitudes miss a substantial fraction of the flux of de Vaucouleurs 
profile galaxies; the Sersic LF is substantially higher than the Petrosian LF at the bright end. 

Key words: galaxies: evolution — galaxies: luminosity function, mass function — galaxies: 
statistics. 


1 INTRODUCTION 


The luminosity function (LF) is perhaps the most fundamental 
model-independent quantity that can be measured from a galaxy 
redshift survey. Reproducing the observed LF is the first require¬ 
ment of a successful model of galaxy formation, and thus accurate 
measurements of the LF are important in constraining the physics 
of galaxy formation and evolution (e.g. [Benson et al.|2003^ . In ad¬ 
dition, accurate knowledge of the survey selection function (and 
hence LF) is required in order to determine the clustering of a flux- 
limited sample of galaxie s (|Cole|2011[). 

A standard 1/Vmax ^Schmidt] 1968^ estimate of the LF is vul¬ 
nerable to radial density variations within the sample. This vul¬ 
nerability can be largely mitigated by multiplying the maximum 
volume in which each galaxy is visible, VLax, by the integrated 


radial overdensity of a density-defining population i Baldry et al.j 
2012|l. Maximum-likelihood methods ( jSandage et al.|1979[ 


2006 


Efstathiou et al.|1988[l, which assume that the luminosity and spa- 
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tial dependence of the galaxy number density are separable, are, 
by construction, insensitive to density fluctuations. However, if the 
sample covers a significant redshift range, galaxy properties (such 
as luminosity) and number density are subject to systematic evolu¬ 
tion with lookback time. All of the above methods must then either 
be applied to restricted redshift subsets of the data, or be modified 
to explicitly allow for evolution (e.g.[Lin et al.|19991|Loveday et al.j 

[20T2I 1. 

jColej ( |201 1[ | recently introduced a joint stepwise maximum 
likelihood (JSWML) method, which jointly fits non-parametric es¬ 
timates of the LF and the galaxy overdensity in radial bins, along 
with an evolution model. In this paper we describe modifications 
made to the JSWML method in order to successfully apply it to the 
Galaxy and Mass Assembly (GAMA) survey ( [Driver et al.|2011) . 
In the GAMA-II sample, L* galaxies can be seen out to redshift 
2 : « 0.35, and so one has a reasonable redshift baseline over which 
to constrain luminosity and density evolution. |Loveday et al.[p012[ ( 
have previously investigated LF evolution in the GAMA-I sample, 
finding that at higher redshifts: all galaxy types were more lumi¬ 
nous, blue galaxies had a higher comoving number density and red 
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galaxies had a lower comoving number density. Here we exploit 
the greater depth (0.4 mag) of GAMA-II versus GAMA-I, and use 
an estimator of galaxy evolution that does not assume a parametric 
form (e.g. a Schechter function) for the LF. 

The paper is organized as follows. In Section]^ we describe 
the GAMA data used along with corrections made for its small 
level of incompleteness. Our adopted evolution model is described 
in Section 1^ and the density-corrected Vinax method in Section]^ 
Methods for determining the evolution parameters are discussed in 
Section!^ We present tests of our methods using simulated data in 
Section]^ and apply them to GAMA data in Section]^ We briefly 
discuss our findings in Sectionj^and conclude in Section]^ 

Throughout, we assume a Hubble constant of Hq = 
lOO/i km s“^ Mpc“^ and an Q,m = 0.3, Ha = 0.7 cosmology 
in calculating distances, co-moving volumes and luminosities. 


2 GAMA-II DATA, K- AND COMPLETENESS 
CORRECTIONS 

In April 2013 the GAMA survey completed spectroscopic cover¬ 
age of the three equatorial fields G09, G12 and G15. In GAMA-II, 
these fields were extended in area to cover 12 x 5 degrees eaclQ 
and all galaxies were targeted to a Galactic-extinction-corrected 
SDSS DR7 Petrosian r-band magnitude limit of r = 19.8 mag. In 
our an^sis, we include all main-survey targets (SURVEYJLASS 
>= djjwith reliable A UTOZ jB aldry et al.|2014 i redshifts {nQ ^ 
3) from TilingCatv43 ( [Baldry e t al.||201 0|. Redshifts (from Dis- 
tancesFramesvl2) are corrected for local flow using the|Tonry et al.| 
( |2000 ^ attractor model as described by|Baldry e t al.| ( [^12^ . 

We calculate LFs using both |Petrosian ^1976^ md |Sersic| 
(|1963| photometry, corrected for Galactic extinction using the dust 
maps of |Schlegel et al.||l998| l. We use single Sersic model magni¬ 
tudes truncated at ten effective radii as fit by |Kelvin et al.| ( [2012| l. 
Kelvin et al. show that these recover essentially all of the flux for 
an n = 1 (exponential) profile, and about 96 per cent of the flux 
of an n = 4 (de Vaucouleurs) profile. SDSS Petrosian magnitudes, 
while also measuring almost all of the flux for exponential profiles, 
measure only about 82 per cent of the flux for de Vaucouleurs pro¬ 
files ( [Blanton et al.|200T| >. Sersic magnitudes are, however, more 
susceptible to contamination from nearby bright objects, which 
can cause them to be overestimated by several mag. We identify 
galaxies which may have contaminated photometry by searching 
for brighter stellar neighbours within a distance, up to a maximum 
of five arcmin, of twice the star’s isophotal radius (isoA_r in 
the SDSS PhotoOb j table). Five per cent of GAMA targets are 
flagged in this way. 

Fig. [T] shows a histogram of Am = rpetro — fSersic for all 
GAMA-II main-survey targets (continuous blue histogram), and 
for targets without a nearby bright stellar neighbour (black dotted 
histogram). The majority (about 72 per cent) of excluded galax¬ 
ies have positive Am, i.e. are brighter in Sersic than Petrosian 
magnitude. The dashed red histogram indicates targets without a 
nearby bright stellar neighbour that are classified as red (as defined 
towards the end of this section). It is clear from this figure that un¬ 
contaminated red galaxies preferentially have brighter Sersic than 


^ The RA, dec ranges of the three fields, all in degrees, are G09: 129.0- 
141.0, -2.0—1-3.0; G12: 174.0-186.0, -3.0—1-2.0; G15; 211.5-223.5, 
-2.0—1-3.0. 

^ Note that in this latest version of TilingCat, objects that failed visual 
inspection (VIS_CLASS = 2,3 or 4) also have SURVEY_CLASS set to zero. 



rprtro-rscrsic (mag) 


Figure 1. Histogram of the difference between Petrosian and Sersic mag¬ 
nitudes for all GAMA-II main-survey targets (continuous blue histogram), 
and for targets without a nearby bright stellar neighbour, as defined in the 
text (black dotted histogram). The dashed red histogram indicates the sub¬ 
set of the latter targets classified as red. The vertical dotted lines denote the 
additional constraint |rpetro — rsersici <2.0 mag required for galaxies to 
be assumed uncontaminated; only about 0.3 per cent of remaining targets 
lie beyond these limits. 


Table 1. Classification of the 487 GAMA targets without a bright stellar 
neighbour (as defined in the text) for which Am = rpetro ~ ’’Sersic lies 
outside the range [—2, 2] mag. See text for meaning of first column. 


Class 

Am > 2 

Am < —2 

OK 

12 

42 

Deblend 

19 

56 

FSC 

13 

22 

BSC 

144 

16 

Merger 

56 

37 

Sky 

17 

46 

NO 

7 

0 

Total 

268 

219 


Petrosian magnitudes. This is as expected, assuming that they are 
bulge-dominated, and hence have profiles with higher Sersic index. 

We exclude an additional 487 targets (0.3 per cent of the to¬ 
tal) for which the r-band Sersic and Petrosian magnitudes differ by 
more than 2 mag. This magnitude difference cut is somewhat arbi¬ 
trary, but is designed to exclude galaxies with bright stellar neigh¬ 
bours that do not quite satisfy the above criterion (for instance if a 
galaxy lies on a star’s diffraction spike) or with bad sky background 
determination. It seems extremely unlikely that the Sersic magni¬ 
tude would recover more flux than this from an uncontaminated 
galaxy. 

We have visually inspected these additional culled targets, for 
which jAm| = |rpetro — J'SersicI > 2 mag, and placed them in 
one of the following categories: OK: no obvious problem; Deblend: 
large galaxy image likely to have been shredded by the SDSS de¬ 
blending algorithm; FSC: nearby faint stellar companion (compara¬ 
ble to or fainter than target); BSC: nearby bright stellar companion 
(much brighter than target); Merger: nearby galaxy companion(s); 
Sky: bad sky background; NO: no object visible. The number of tar¬ 
gets falling into each category, subdivided by whether Am is pos- 
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itive (Sersic flux is brighter) or negative (Petrosian flux is brighter) 
is given in Table[^ For the former sample, just over half of the cases 
of possibly overestimated Sersic flux appear to be due to a nearby 
star which has more successfully been excluded from the Petrosian 
flux estimate. For the latter sample, the most common cause of un¬ 
derestimated Sersic flux or overestimated Petrosian flux is likely 
due to deblending issues or a bad sky determination. We note that 
the presence of a nearby bright star should be totally uncorrelated 
with a galaxy’s intrinsic properties, and so excluding targets for this 
reason should not bias the sample in any way. A small bias could 
be caused by excluding the ~ 20 per cent of inspected galaxies 
(about 0.06 per cent of total targets) for which the suspect photom¬ 
etry is caused by a neighbouring galaxy, since galaxies in crowded 
regions are expected to be more luminous than average. In cases 
when the Petrosian and Sersic magnitudes differ by more than 2 
mag, both magnitude estimates are suspect, and so it is debatable 
whether these objects should be GAMA targets at all. At worst, 
the effect of excluding targets with bright stellar neighbours or dis¬ 
crepant magnitudes (5.3 per cent of the entire GAMA-II sample) 
will be to bias the LF normalization low by up to five per cent. 

After excluding GAMA main survey targets with either an un¬ 
reliable redshift (1.2 per cent) or suspect photometry (5.3 per cent), 
we are left with a sample of 173,527 galaxies in the redshift range 
0.002 < « < 0.65. 


To determine if -corrections, we use KCORRECT v4.2 an-| 
[ton & Roweis| |2007|) to fit spectral energy distributions to ugriz 
GAMA matched-aperture SExtractor jBertin & Amouts [1996^ 
AUTO magnitudes taken from ApMatchedCatv04 ( [Hill et al.| 
|2011|l. As shown in Appendix B of [Taylor et al.H20lT) and Fig. 17 
of Kelvin et al.j j2012|l, SDSS model magnitudes, which have 
been recommended for calculating galaxy colours, e.g. |Stoughtolt| 
jet al.j {2002) , are ill-behaved for galaxies of intermediate Sersic 
index which are well fit by neither pure exponential nor pure de 
Vaucouleurs profiles. GAMA matched-aperture magnitudes do not 
force a particular functional form on the galaxy profile and so pro¬ 
vide more reliable colours for all galaxy types. In practice, we find 
that the choice of magnitude type used for TT-corrections makes lit¬ 
tle difference to our LF estimates, with the Schechter fit parameters 
changing by less than 1-sigma. We use iT-corrections to reference 
redshift zo = 0.1 in order to allow direct comparison with previous 
results ( jLoveday et al.|20r^ . For the three GAMA-II targets which 
are missing AUTO magnitudes, and for the 3.3 per cent of targets 
for which KCORRECT reports a statistic of lO.O or larger, im¬ 
plying a poor SED fit, we set the TT-correction to the mean of the 
remaining sample. We have visually inspected 235 of these targets 
with poor-fitting SEDs. About 29 per cent are close to a bright star 
or are otherwise likely to suffer from poorly-estimated sky back¬ 
ground; about 22 per cent have one or more close neighbours and 
may thus suffer contaminated photometry; about 13 per cent show 
evidence of AGN activity. The remaining 35 per cent show no ob¬ 
vious reason for the SED fit to be poor, but it seems likely that 
many of these cases may be due to poor u-band photometry with 
underestimated errors. 

While SDSS DR7 has improved photometric calibration over 
DR6 (used for selection of GAMA-I targets), it will suffer from the 
same surface-brightness-dependent selection effects as DR6, and so 
we assume the same imaging completeness Cim as shown in Eig. 1 
of jLoveday et al.|(2012| >. In this paper, we only measure the r-band 
LE, and so assume that target completeness is 100 per cent. In fact, 
just 0.1 per cent of GAMA-II main targets with r < 19.8 mag lack 
a measured spectrum, with no systematic dependence on magnitude 
(Liske et al., submitted to MNRAS). Since GAMA-II uses a new. 




Figure 2. Redshift success rate as a function of r-band fibre magnitude. The 
top panel shows histograms of rgbre for all observed galaxies in blue, and 
galaxies with a reliable redshift measurement (nQ > 2) in green. Redshift 
success, the ratio of the latter to the former, is shown as a histogram in the 
lower panel, along with a best-fit sigmoid-type function. The large fluctu¬ 
ations at faint magnitudes (rgure > 21) are simply due to small-number 
statistics: the success rate is the ratio of two small numbers. 



log^o weight 


Figure 3. Histogram of completeness-correction weights for GAMA-II 
galaxies. Note that both axes use logarithmic binning. The vertical line cor¬ 
responds to the applied upper limit weight cap of 5.0. 


fully-automated redshift measurement ( [Baldry et al.|2014) , we have 
re-assessed redshift success rate for GAMA-II. Fig. shows red¬ 
shift success rate, defined as the fraction of observed galaxies with 
reliable (nQ ^ 3) redshifts, as a function of r-band fibre magni¬ 
tude. This success rate is well-fit by a modified sigmoid function 

a = (1) 

with parameters a = 2.55 mag“^, b = 22.42 mag and c = 2.24. 
The extra parameter c (c.f. [Ellis & Bland-Flawthom|2007||Loveday| 
jet al.|201^ is introduced to provide a more extended decline in Cz 
around rnbre ~ 20 mag. Without it, the sigmoid function drops too 
sharply to faithfully follow the observed Cz ■ 

Each galaxy is given a weight equal to the reciprocal of the 
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product of imaging completeness and redshift success rate, Wi = 
l/(C'imiC'zi). A histogram of these weights is shown in Fig. 
While the vast majority of galaxies (99.5 per cent) have Wi < 2, 
there is a tail of rare objects with weights as high as 100 or more. 
We have visually inspected the 157 objects with an assigned weight 
above 10.0. Of these, 38 per cent are close to a bright star or are oth¬ 
erwise likely to have a poorly-determined sky background; another 
38 per cent have nearby neighbouring galaxies which might lead to 
a compromised surface-brightness estimate; 10 per cent are isolated 
and show no obvious visual indication of being of low surface- 
brightness. That left just 14 per cent which appeared to be genuine 
low surface-brightness galaxies, potentially with half-light surface 
brightness /rso.r ^ 24 mag arcsec“^ and/or with fibre magnitude 
fflbre ^ 22 mag. We therefore chose to set an upper limit cap of 5.0 
on incompleteness weights, i.e. to set Wi — min(Wi,5.0). This 
limit corresponds to the inverse redshift success rate for galaxies 
with the faintest fibre magnitudes (Fig.|^. While only 297 galaxies 
(0.16 per cent of the total) have Wi > 5.0, these galaxies are likely 
to he at the extreme faint end of the LF, where there are few ob¬ 
served galaxies, and so spurious weights could potentially bias the 
LF faint end. The mean galaxy weights before and after applying 
this cap are 1.12 and 1.09, respectively. 

The effect of applying this weight cap is to reduce the best-fit 
value of the density evolution parameter Pe by about 40 per cent, 
with a corresponding increase in the best-fit value of the luminosity 
evolution parameter Qe- Best-fit LF parameters change by less than 

1-(T. 

When subdividing GAMA galaxies into blue and red sub¬ 
samples, we use the colour cut of [Loveday et al.| ( |2012| ), namely 

“■^(5-f)Kron =0.15-0.03°'^(M, -51og/l). (2) 

A detailed investigation of colour bimodality in GAMA has re¬ 
cently been presented by |Taylor et al.|p015^ . They utilise restframe 
and dust-corrected [g — i) colour, and argue that a probablistic as¬ 
signment of galaxies to ’R’ and ’B’ populations is preferable to 
a hard (and somewhat arbitrary) red/blue cut. They also empha¬ 
sise that colour is not synonymous with morphological type, but 
rather provides a proxy for mean stellar age within a galaxy. Also, 
of course, a galaxy may appear red in uncorrected restframe colour 
due to dust extinction, rather than an old stellar population. In this 
paper, we stick with the simple colour-cut of equation for two 
reasons: (i) to allow direct comparison with the results of[Loveday| 
|et al.| ( [20T^ ; (ii) the |Taylor et al.| ( [2OT3| > model of the colour-mass 
distribution has been tuned to a nearly volume-limited sample of 
galaxies at redshift z < 0.12 — the model parameters are likely to 
evolve at higher redshift. 

Uncertainties in measured quantities, such as radial overden¬ 
sity and the LF, are determined by jackknife resampling. We sub¬ 
divide the GAMA-II area into nine 4x5 degree regions, and then 
recalculate the quantity nine times, omitting each region in turn. 
For any quantity x, we may then determine its variance using 

\ax{x)=^^-j^^{xi-xf, (3) 

i=l 

where A = 9 is the number of jackknife regions, Xi is our estimate 
of X obtained when omitting region i, and x is the mean of the 
Xi. The numerator {N — 1) in the pre-factor allows for the fact 
that the jackknife estimates are not independent. Each jackknife 
region contains an average of 19,281 galaxies for the full GAMA- 
II sample (i.e. without colour selection). 


3 PARAMETRIZING THE EVOLUTION 

We parametrize luminosity and density evolution over the redshift 
range 0.002 < z < 0.65 using the parameters Qe and Pe in¬ 
troduced by |Lin et ah] ( |1999| l. This model assumes that galaxy 
populations evolve linearly with redshift in absolute magnitude, 
parametrized by Q^, and in log number density, parametrized by 
Pe. Specifically, the luminosity e-correction is given by E(z) — 
Qe{z — zo), such that absolute magnitude M is determined from 
apparent magnitude m using 

M = m- 51ogi(, di(z) - 25 - K{z-, zq) -F Qe{z - zq), (4) 

where dL{z) is the luminosity distance (assuming the cosmolog¬ 
ical parameters specified in the Introduction) at redshift z and 
K{z; zo) is the P-correction, relative to a passband blueshifted by 
Zq. Luminosity evolution is determined relative to the same redshift 
zo = 0.1 as the A-correction. 

Evolution in number density P{z) is parametrized as 

P{z) = = P{z = 0)10“'‘^=L (5) 

The motivation for this choice of parametrization is that if the shape 
of the LF does not evolve with redshift, that is it shifts only hori¬ 
zontally in absolute magnitude by Qe, and vertically in log-density 
by Pe, then luminosity density px evolves as 

Pl{z) = (6) 

While Pe and Qe are strongly degenerate, and so poorly con¬ 
strained individually, their sum Pe + Qe is well-constrained 
|et al. 1 199^ [Loveday et al.|2012| l. We set further constraints on the 
linear combination of these parameters in Section[7] 


4 DENSITY-CORRECTED Umax METHOD 

In this section we describe our technique for determining the LF 
using a maximum-likelihood, density-corrected Umax estimator, as¬ 
suming that evolution is known. We will discuss how we determine 
the evolution parameters Qe and Pe in Section Our method is 
based on the joint stepwise maximum likelihood (JSWML) method 
of [Cole[ ( [201 l[ l, which jointly fits the LF and overdensities in radial 
bins of redshift caused by large-scale structure. Cole’s derivation 
starts with an expression for the joint probability of finding a galaxy 
at specified redshift and luminosity, and assumes that all galaxies 
have identical evolution- and iT-corrections. We wish to allow for 
individual K- (and in the future e-) corrections, in which case it 
is easier to start with the conditional probability that an observed 
galaxy of luminosity Li has a redshift Zi, assuming that the lu¬ 
minosity and spatial dependence of the galaxy number density are 
separable. This conditional probability is given by ([Saunders et al.j 

|T^ : 

F* A(z)P(z)^dz' 

Here we have factored the mean density at redshift z, n(z) — 
A(z)P(z), into a product of the galaxy overdensit)]^ A(z) due to 
large-scale structure times the steadily evolving density P(z) from 
equation ij^; dVjdz is the differential of the survey volume, and 
2max,i IS the maximum redshift at which galaxy i would still be 

^ Following [Cole[ f20 11 we use the term overdensity to mean a multi¬ 
plicative relative density, so that A = 1 corresponds to average density. 
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visible, determined by the survey flux limit along with the galaxy’s 
luminosity, K- and e-corrections. 

Adopting binned estimates of the galaxy overdensity A, and 
weighting each galaxy by its incompleteness-correction weight, 
Wi, we obtain a log-likelihood 


the python module pYqt_fit. kd^ does not, of course, 
correspond to the true galaxy LF, but rather to the LF convolved 
with a Gaussian of standard deviation b. Therefore when plotting 
the LF and fitting a Schechter function, we use the standard binned 
LF rather than 


ln£ = ^ VFi 

i 


In y ^ Aj Pj Vj Dij 
j 


3 


( 8 ) 


5 DETERMINING EVOLUTION PARAMETERS 


Flere Vj , Pj and Aj are the volume, density evolution and galaxy 
overdensity respectively in redshift bin j ; the function Dij is a sim¬ 
ple binning function, equal to unity if galaxy i lies in redshift bin 
j, zero otherwise, and Sij is the fraction of redshift bin j in which 
galaxy i is visible. In the present analysis we employ redshift bins 
of width A 2 = 0.01. The maximum-likelihood solution for the 
overdensities Aj, given by dln£/dAj — 0, may be obtained by 
iteration from: 


Aj — IVsum,j 


E 


y m 


(9) 


where VFsumj = WiDij is the sum of galaxy weights in red¬ 
shift bin j and UmL* = AkPkVkSik, the effective volume, 
corrected for evolution and fluctuations in radial density, within 
which galaxy i is visible. 

The LF, unaffected by density fluctuations, may then be esti¬ 
mated by substituting the usual expression for VmsLX- 


= E 


W^Dg 


( 10 ) 


where Dg = 1 if galaxy i is in luminosity bin I, zero other¬ 
wise. |Cole| pOl 1[ | shows that this expression may be derived via 
maximum likelihood, at least in the case of identical e- and K- 
corrections. 

Cole also discusses an extension to this method whereby pa- 
rameter(s) describing the density evolution P{z) may be deter¬ 
mined simultaneously with the overdensities Aj by adding prior 
constraints on the values of Aj using the known clustering of galax¬ 
ies. However, for our choice of density evolution parametrization 
(equation |^, the derivative in Cole equation (25) no longer de¬ 
pends explicitly on the evolution parameter, leading to a lack of 
convergence. We therefore prefer to search over both luminosity 
and density evolution parameters, as described in the next section. 

A stepwise estimate of the LF, as given by equation l |10[ l, is not 
constrained to vary smoothly from bin to bin. Furthermore, at very 
low and high luminosity there may be bins containing no galaxies, 
resulting in an ill-defined log-likelihood (see equation [T^ below). 
This problem is exacerbated when exploring possible values of the 
luminosity evolution parameter Qe, as galaxies will then shift from 
bin to bin as Qe is varied, resulting in unphysical sharp jumps in 
likelihood. To overcome these problems, we employ a Gaussian- 
smoothed estimate of the LF: 



Here the smoothing kernel G is a standard Gaussian, b is the 
smoothing bandwidth. Mi is the (K- and e-corrected) absolute 
magnitude of galaxy i and Mi is the absolute magnitude at the cen¬ 
tre of bin 1. In order not to underestimate the extreme faint-end of 
the LF, it is important to apply boundary conditions to corre¬ 
sponding to the chosen range of absolute magnitudes. We do this 
using the default renormalization method and bandwidth choice of 


In Cole’s original derivation of this method, one maximises a pos¬ 
terior likelihood (Cole equation SsQover the luminosity evolution 
parameter Qe (Cole calls this parameter u). When applying this 
method to GAMA data, we found that the estimated value of Qe 
diverged, unless one places an extremely tight prior on its valu^ 
Our problem was traced to the fact that varying Qe changes all 
of the inferred absolute magnitudes (as well as visibility limits) 
for each galaxy. Choosing fixed absolute magnitude limits within 
which to determine the LF thus results in a change of sample size 
as Qe varies, leading to likelihoods that cannot be directly com¬ 
pared. Even if one includes the term on the second line of Cole 
equation (36), which yields — Wtot In Atot in the case of identical 
K- and e-corrections, the estimate of Qe still diverges as galaxies 
shift systematically brighter or fainter as Qe decreases or increases. 
We therefore consider two alternative methods to optimize the evo¬ 
lution parameters. 


5.1 Mean probability 

Our first solution is to consider not the product of the probabilities 
of observing each galaxy, but instead the geometric mean of the 
probabilities which does not vary systematically with sample size 
N. Our pseudo-log-likelihood In V is then given by (Cole equation 
36) 

E hvjPAj) + ^ E 

3 ^ 

- ^ E E E |L0 

i j I 

^ [Aj - If {Qe - Qof jPe - Pof 

j 20-^^ 2(Tq^ 2(Tp^ 

(12) 

where VFaum.j is the sum of galaxy weights in redshift bin j, 
H4um,i is the sum of galaxy weights in luminosity bin I, and 
S{Luiin,i,j\Li) is the fraction of luminosity bin I for which galaxy i 
at redshift Zj would be visible. The term on the second line is a con¬ 
stant in the case of identical K- and e-corrections; with identical 
e- but independent A'-corrections we find that including this term 
makes a negligible difference to the maximum-likelihood solution. 


https : //pypi . python . org/pypi/PyQt-Fit 
® Note that Cole equations (36-38) are missing factors of Pp, such that 

each occurrence of Vp should read VpPp. _ 

® We believe that the reason that the test described in Section 5 of IColel 
(MT) was successful was due to (i) placing a very tight prior {a = 0.05) 
on the density evolution parameter, and (ii) simulating a very deep galaxy 
survey (extending to magnitude r < 24 and redshift 2 : < 1.4). Both of 
these factors minimize the degeneracy between luminosity and density evo¬ 
lution, and hence aid convergence. The GAMA-II sample is signihcantly 
shallower (r < 19.8, z < 0.65), and we do not wish to place tight prior 
constraints on either of the evolution parameters. 
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The terms on the third line are priors on the radial overdensities 
Aj and the evolution parameters Pe and Qe - The priors on are 
essential, as these values are completely degenerate with the den¬ 
sity evolution parameter P,- As discussed by Cole, the expected 
variance in Aj is given by 




1 -I- A-Khj J 3 
hjVj 


(13) 


with hj the predicted density and Vj the volume of redshift bin 
j. The factor J 3 = J r^^{r)dr « 2 , 000 h,“®Mpc® accounts for 
the fact that because galaxies are clustered, they tend to come in 
clumps of 47 rn J 3 galaxies at a time ^Peebles 1980| ). We find, how¬ 
ever, that much more reliable estimates of cr a . are obtained from 
jackknife sampling — see Fig. ^ below. This is particularly true 
in the higher redshift bins, where one is sampling the clustering of 
the most luminous galaxies, and where adopting a universal value 
for J 3 underestimates the actual density fluctuations observed be¬ 
tween jackknife samples. The priors on Pe and Qe are optional, 
and may help convergence in some cases. We adopt broad pri¬ 
ors of {Qo,^q^) = (1,1) and (Poj'^p^) ~ (2,1). These values 
were chosen to be consistent with the findings of [Loveday et al.| 
( | 2012 [ while still allowing some freedom for the optimum values 
to change under the present analysis. 


5.2 LF-redshift 


Our second method compares LFs estimated in two or more red¬ 
shift ranges: if the evolution and density variations are correctly 
modelled, then the LFs should be in good agreement; if evolution 
parameters are poorly estimated, then one would expect poor agree¬ 
ment. We then minimize the x^ (= — 2 In £) given by 


X 


2 


I ) + Var((/.f) Y 


(14) 


where 0; is the Gaussian-smoothed LF in magnitude bin I for the 
broad redshift range j, and Var((^j) is the corresponding variance, 
determined by jackknife resampling. We restrict the sum over mag¬ 
nitude bins I to those bins which are complete given the redshift 
limits (see Section 3.3 of jLoveday et al.|2012^ and which include 
at least ten galaxies for all values of Qe between specified lim¬ 
its. In practice, we have found best results are achieved using just 
two redshift ranges, split near the median redshift of the sample, 
a « 0.2, so that the ‘knee’ region of the LF around L* is well- 
sampled by both, and hence the degeneracies between luminosity 
and density evolution are minimized. If one chooses three or more 
redshift ranges, there will be very little luminosity coverage in com¬ 
mon to the lowest and highest ranges, and so one does not really 
gain much information in doing so. Again, it is essential to place a 
prior on the overdensities (final sum in equation 14 with cr^ also 


determined from jackknife resampling) to remove the degeneracy 
with density evolution. This method places no priors on the values 
for the evolution parameters. 


5.3 Finding optimum evolution parameters 

We first evaluate x^ values, using each of the above methods, on 
a rectangular grid of (Pe, Qe), thus allowing one to visualise the 
correlations between the evolution parameters. The grid point with 
the smallest x^ value is then used as a starting point for a downhill 
simplex minimisation to refine the parameter values corresponding 
to minimum x^- 


In order to quantify the degeneracy between evolution param¬ 
eters, we slice the x^ grid in bins of Pe. For each slice we fit a 
quadratic function to X^iQe) using the five {Qe, X^) values clos¬ 
est to the point of minimum x^ in that slice. Using this quadratic 
fit, we locate the point ^2 of minimum x^ and its 1 -sigma 
range, i.e. the range of Qe values where x^ increases by unity from 
the minimum. We find both for simulations and for real data that 
the Qj, ^2 -Pe relation is very well fit by a straight line, and so 
we perform a linear least-squares fit to (Pe, Qe ^2 ) to obtain the 
relation Qe = mPe -h c which minimizes x^- 

6 TESTS USING SIMULATED DATA 
6.1 The simulations 

In this section, we test our implementation of the IS WML estimator 
using simulated data, following the procedure outlined in Section 5 
of [CoIe| ( [20TT] l. 

We start by choosing a model LF with jSchechterj and 

evolution parameters close to those obtained from the GAMA-I sur¬ 
vey by [Loveday et al.j ( [2012) and as given in Table We then 
randomly generate redshifts with a uniform density in comoving 
coordinates, modulated by our assumed density evolution (equa¬ 
tion]^, over the range 0.002 < 2 < 0.65. Absolute magnitudes 
are selected randomly according to our assumed Schechter func¬ 
tion from the range —24 < M < —12. From each absolute mag¬ 
nitude we subtract Qe{z — 0.1) to model luminosity evolution. We 
then assign apparent magnitudes r using Ff-correction coefficients 
selected randomly from the GAMA-II data and reject simulated 
galaxies fainter than r = 19.8. This process is repeated until suf¬ 
ficient random galaxies have been generated to give the required 
number density, 

r^max r r.max i^max) WU 

/ / 4>{L,z)dL — dz, (15) 

within a volume corresponding to that of the three GAMA-II fields, 
viz 3 X 5° X 12° = 180 deg^ 

In order to simulate the effects of galaxy clustering, we spilt 
the simulated volume into 65 redshift shells p of equal thickness 
A 2 « 0.01 and with volume Vj. For each shell we generate a ran¬ 
dom density perturbation Sj drawn from a Gaussian with zero mean 
and variance AnJalVj, with TrrJa = 30, 000/i“^Mpc®. We then 
randomly resample A)) = (1 -|- Sj)Nj of the original Nj simulated 
galaxies in each shell p, thus producing fluctuations consistent with 
the assumed value of J 3 . 

Imaging completeness and redshift success are modelled by 
generating surface brightnesses and fibre magnitude for each sim¬ 
ulated galaxy according to the relations observed in GAMA-I data, 
see Appendix Al of [Loveday et al.[ ( [2012). Imaging co mpleteness 
Cini is then determined from Fig. 1 of [Loveday et ah] ( [2012) and 
redshift success Cz from equation Q. Simulated galaxies are then 
chosen randomly with probability equal to CimCz and assigned 
a weight Wi = l/(CimiC'zi) to compensate for those simulated 
galaxies omitted from the sample. 

This procedure is repeated to generate ten independent mock 
catalogues, each containing around 180,000 galaxies. These mock 
catalogues are run through the JSWML estimator, with evolution 
parameters being determined using both methods discussed in the 
previous section. Since the mock galaxies are clustered only in red¬ 
shift shells and not in projected coordinates on the sky, we de¬ 
termine the expected variance in overdensity using equation 
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Table 2. Mean and standard deviation of the evolution and Schechter pa¬ 
rameters recovered from ten simulated GAMA catalogues. Parameters m 
and c quantify the linear relation Qe = mPe -P c which minimizes . 



True 

mean prob 

LF-redshift 

Qe 

0.7 

0.68 ± 0.23 

0.61 ± 0.14 

Pe 

1.8 

1.71 ± 0.67 

1.88 ±0.54 

Cov{Qe,Pe) 


-0.16 

-0.08 

m 


-0.37 ±0.01 

-0.33 ±0.01 

c 


1.31 ± 0.09 

1.24 ±0.07 

a 

-1.23 

-1.23 ±0.02 

-1.23 ±0.01 

M* - 5 Ig /t 

-20.70 

-20.72 ± 0.06 

-20.74 ±0.03 

lg((/i*/^^Mpc-3) 

-2.00 

-2.10 ±0.09 

-2.12 ±0.07 



Figure 4. 95 per cent confidence limits on evolution parameters determined 
from ten simulated datasets (light contours) and their average (heavy con¬ 
tour) determined using (t op) mean probability (equation 1 12) and (bottom) 
LF-redshift (equation |14| . The error bars show the mean and standard 
deviation of the (Pe,(3e) parameters from each simulation which yield 
minimum x^- The input evolution parameters for these simulations were 
Pe = 1.8, Qe = 0.7. 

rather than jackknife resampling. We employ 65 redshift shells out 
to z = 0.65 and calculate the LF in bins of AM — 0.25 mag over 
the range —23 < M < —15 mag. 

6.2 Simulation results 

The mean and standard deviation of each recovered parameter, and 
the covariance between evolution parameters, are given in Table|^ 
We see that the input evolution and LF parameters are recovered 
within about one standard deviation for both methods. 

Fig. El shows 95 per cent confidence limits on the evolution 
parameters measured from each of the simulations. We see that 


the error contours are significantly smaller using the LF-redshift 
method compared with the mean probability method. Flowever, 
this test is idealized, in that our choice of evolution parametriza- 
tion is identical in the simulations and in the analysiQ and so we 
will apply both methods of constraining evolution parameters to 
the GAMA data in the following Section. Note that the simulations 
have no inbuilt covariance between evolution parameters: they all 
use identical values of Pe and Qe- The degeneracies (as quantified 
by Cov(Qe, Pe) and the parameters m and c in Table|^ arise as 
a result of the fitting process. For an LF described by an unbroken 
power law, the degeneracy between Pe and Qe would be total, i.e. 
evolution in luminosity and density would be indistinguishable. 

7 RESULTS FROM GAMA 

7.1 Evolution 

Fig.|^shows 95 per cent confidence limits on the evolution param¬ 
eters Pe, Qe determined using equations l |12| l and ( |14) for the full 
GAMA-II sample and for blue and red galaxies separately. We see 
that the confidence limits obtained with the two different methods 
largely overlap, although there are small differences between them. 
Best fit evolution parameters are given in Table The difference 
in LFs obtained using evolution parameters determined with the 
two different methods is negligible (much less than the 1 -a random 
errors; see Table |^. This illustrates the robustness of the LF es¬ 
timate to the individual values assumed for P^ and Qei as long as 
their joint estimate is reasonable, e.g. they lie within the 95 per cent 
likelihood contours of Fig.then overestimating one evolution pa¬ 
rameter (e.g. Pe) is largely compensated for by underestimating the 
other (e.g. Qe). 

The differences in density evolution (Pe) for red and blue 
galaxies are not significant. Blue galaxies do however exhibit sig¬ 
nificantly stronger evolution in luminosity (Qe) and in luminosity 
density (Qe + Pe) than red galaxies, at the ~ 5-a level. 

The differences between red and blue galaxies agree qualita¬ 
tively with those of |Loveday et al.| ( |2012) , although in the present 
analysis we no longer see any evidence for negative density evolu¬ 
tion for red galaxies. The three samples show very similar degen¬ 
eracies in (Pe, Qe) parameter space. The errors on Pe and Qe in 
Table are the formal errors obtained by holding one parameter 
fixed and varying the other until increases by one. Given the 
scatter in 95 per cent confidence limits between simulations shown 
in Fig. 12 more realistic errors, and their covariance, may be ob¬ 
tained from Table 1^ 

Since the exact values assumed for the evolution parameters 
have such a small effect on the LF parameters, see Tablej^below, 
for the remainder of this paper we assume evolution parameters 
found from the LF-redshift method in the lower half of Table[3] 

7.2 Radial overdensities 

Radial overdensities are shown in Fig. |2 While our evolution 
model is performing well, insofar as A( 2 ) oscillates about unity, for 
redshifts 2 : < 0.5, beyond this limit the overdensities are systemat¬ 
ically high. This effect is almost entirely due to red galaxies, sug¬ 
gesting that luminosity and/or density evolution increases sharply 
at 2 « 0.5 for these galaxies compared with our model (sectionj^l. 

^ We are performing a self-consistency test. It is unlikely that real galaxy 
populations evolve exactly according to our parametrization. 


© 0000 RAS, MNRAS 000, 000-000 













8 J. Loveday et al. 


Table 3. Best-fitting evolution parameters for GAMA-II galaxy samples obtained using both mean probability and LF-r edshift methods. Parameters m and c 
quantify the linear relation Qe = mPe + c. For the LF-redshift method only, is the reduced x^ from equation (14); The uncertainties quoted on Qe and 
Pe come from the bounding box containing the I-sigma likelihood contour; the uncertainty on Pe + Qe is given by the distance from the point of minimum 
to the 1-sigma likelihood contour along the direction Pe = Qe- 


Sample 

Qe 

Pe 

Qe Pe 

m 

c 

xi 

Mean Probability 

All 

1.03 ±0.10 

1.00 ±0.25 

2.02 ±0.05 

-0.36 

1.38 


Blue 

1.09 ±0.10 

1.30 ±0.25 

2.39 ±0.04 

-0.35 

1.55 


Red 

0.58 ±0.18 

1.55 ±0.40 

2.12 ±0.08 

-0.38 

1.17 


LF-redshift 

All 

1.03 ±0.07 

1.00 ±0.20 

2.02 ±0.05 

-0.35 

1.37 

3.76 

Blue 

1.18 ±0.05 

1.07 ± 0.15 

2.25 ±0.04 

-0.34 

1.55 

3.46 

Red 

0.73 ±0.10 

1.25 ±0.25 

1.98 ±0.06 

-0.36 

1.16 

3.35 



Figure 5. 95 per cent confidence limits on GAMA-II evolution parameters 
for all, blue and red galaxies as labeled. The upper panel shows the limits 
obtained using mean probability (eq uatio n |12) ; the lower panel shows re¬ 
sults using LF-redshift (equation |l4| . The large dots indicate the loca¬ 
tion of minimum x^ ■ The large errorbars show the evolution parameters and 
68 per cent confidence limits estimated for the combined GAMA-I sample 
in the r band by |Loveday et al.H2012| Qpar and Ppar from Table 5). 

It seems unlikely that incompleteness corrections could cause this, 
as there is no noticeable increase in weights beyond a = 0.5. Only 
0.8 per cent of GAMA-II main survey galaxies lie beyond 2 = 0.5, 
too few to constrain a more complicated evolution model, or to look 
for a large overdensity at these redshifts, see Fig.|^ 

Below redshifts 2 = 0.5, we see the same features in radial 
overdensity in all three samples, although the fluctuations, as ex- 


Table 4. Best-fitting r-band LF parameters for GAMA-II galaxy samples 
obtained using both mean probability and LF-redshift methods. For the 
latter method, we show LF parameters obtained using both Petrosian and 
Sersic magnitudes, xi i® the reduced x^ from least-squares Schechter func¬ 
tion fits to the LF estimates; none of the LFs are well-fit in detail by a 
Schechter function, particularly at the bright end. The uncertainties quoted 
on the LF parameters come from jackknife sampling, but do not explicitly 
include the large degeneracies between them. 


Sample 

a 

M* — 5 log h 

log0*/fi^Mpc ^ 

xl 

Mean Probability Petrosian 

All 

-1.26 ±0.07 

-20.71 ± 0.05 

- 2.02 ±0.04 

2.33 

Blue 

-1.38 ±0.06 

-20.36 ±0.05 

-2.27 ± 0.05 

0.96 

Red 

-0.79 ±0.11 

-20.68 ±0.06 

-2.23 ±0.05 

3.48 

LF-redshift Petrosian 

All 

-1.26 ±0.07 

-20.71 ± 0.05 

-2.02 ±0.04 

2.33 

Blue 

-1.37 ± 0.06 

-20.35 ±0.05 

-2.24 ± 0.05 

1.01 

Red 

-0.77 ± 0.11 

-20.64 ± 0.05 

-2.20 ± 0.04 

3.01 

LF-redshift Sersic 

All 

-1.30 ±0.06 

-20.88 ±0.06 

-2.13 ±0.04 

2.96 

Blue 

-1.39 ±0.07 

-20.40 ±0.06 

-2.27 ± 0.05 

1.03 

Red 

-0.79 ±0.12 

-20.72 ±0.07 

-2.21 ± 0.05 

4.91 


pected, are slightly more pronounced in the red galaxy sample. 
Note that the error bands given by equation l |13| > (shaded regions 
in Fig. 1^ are significantly larger/smaller than the jackknife errors 
at low/high redshift. There are two reasons for this: (i) the low- 
redshift bins sample too small a volume for the Js integral to have 
converged, and (ii) the low/high-redshift bins are dominated by 
faint/luminous galaxies, with weaker/stronger clustering than the 
average defined by the assumed value of J3. This is why we use 
jackknife errors rather than the predicted variance in determining 
cr|^. We have tried halving the number of redshift bins to 32, ver¬ 
ifying that the fitted parameters are insensitive to the redshift bin¬ 
ning, with parameters changing by less than one sigma when the 
redshift bin size is doubled from Az = 0.01 to A 2 = 0.02. 

7.3 LFs 

Petrosian and Sersic r-band LFs are shown in Fig. Surface 
brightness and redshift incompleteness have been taken into ac¬ 
count by appropriately weighting each galaxy (up to a maximum 
weight of 5.0, Section)^. We have fit a Schechter function to each 
binned LF using least squares; the fit parameters are tabulated in 
Table 1^ Note that the Schechter fit for red galaxies underestimates 
the faint end of the LF (as well as the bright end — see below). It is 
likely that the faint-end upturn for red galaxies is at least partly due 
to the inclusion of dusty spirals in this sample; the luminosity and 
stellar mass functions of E-Sa galaxies of|Kelvin et al.|p014a|b|l 
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Redshift 2 

Figure 6. Radial overdensities determined from GAMA-II using the entire 
sample and blue and red subsets as labelled, assuming evolution parameters 
as given in the lower half of Table The error bars show uncertainties 
estimated from jackknife sampling and the shaded regions centred on A = 
1 show the expected variance from equation (HI- 


show no indication of a faint-end or low-mass upturn. Fig. 5 of 
IKelvin et al.| ( |2014bt shows that while very few galaxies with el¬ 
liptical morphology are blue ((g — i)o < 0.6), the converse is not 
true: a substantial number of galaxies with spiral morphology are 
red ({g — i)o > 0.8). Any upturn in the luminosity or mass func¬ 
tion of spheroidal galaxies is more likely to be due to the presence 
of so-called little blue spheroids jKelvin et al.|2014b| Fig. Al). Fi¬ 
nally, we note that |Taylor et al.| ( |2015[ ) have shown that the shape 
of the low-mass end of the stellar mass function of red galaxies is 
sensitive to how ’red’ is defined. A low-mass upturn is seen when 
using the definition of|Peng et al.H2010), but not when using those 
of |Bell et al.lpOO^ and |Baldry et al.| ( |2004| l. 

The red galaxy LF, and that for the combined sample, show 
a bright-end excess: there are significantly more high-luminosity 
(Mr — 5 Ig /i < —23 mag) galaxies than predicted by the Schechter 
function fit. This is particularly true for the LF measured using 
Sersic magnitudes, which capture a larger fraction of the total light 
for de Vaucouleurs profile galaxies which dominate the bright end 
of the LF (e.g. |Bemardi et al.|[2013[ l. A bright-end excess above 
a best-fitting Schechter function has been observed in many other 
surveys (e.g. Loveday et al.|199^ [Norberg et al.|2002[ |Montero-| 


|Dorta & Prada||2009^ and appears to be particula rly pronounced 

in bluer bands (e.g. Montero-Dorta & Prada|200q [Loveday et al.| 


© 0000 RAS, MNRAS 000, 000-000 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Redshift 


Figure 7. Redshift histograms for the whole GAMA-II sample and for blue 
and red galaxies separately. The curves in each panel give the predicted 
redshift distribution based on our evolving LF model fits. 


|2012[ [Driver et al.|2013| l. As [Driver et al.| ( |2()13| l point out, given 
the approximately Gaussian distribution of galaxy colours, the LF 
cannot be well fit by a Schechter function in all bands. One should 
however be aware of the possibility that Sersic magnitudes, extrap¬ 
olated as they are out to ten effective radii, are susceptible to over- 
(or under-) estimating the flux of even isolated galaxies if the Sersic 
parameters are poorly fit (although the fitting pipeline does attempt 
to trap for poor fits). Lienee we also show LFs using more stable 
Petrosian magnitudes. 

Our Schechter fits to these LFs are consistent with the r- 
band LFs determined from the GAMA-I sample by [Loveday et al.[ 
( [2012^ , using slightly different methods, and shown in Fig. as 
dotted lines. We also show the ‘corrected’ LF from the IBlantonl 
[et al.[ ( |2005l l low-redshift SDSS sample (without colour selection). 
Considering that this plot is comparing the LFs of SDSS galaxies 
within only 150/i“^Mpc with GAMA galaxies out to 2 « 0.65, the 
agreement is remarkably good, and provides further evidence that 
the simple evolutionary model adopted allows one to accurately 
recover the evolution-corrected LF, despite its poor performance 
beyond redshift 2 ~ 0.5 (Fig.j^. 


7.4 Testing the evolution model 

In Fig.j^ we investigate how faithfully our simple evolution model, 
namely one in which log-luminosity and log-density evolve lin- 





















































































































10 


J. Loveday et al. 



0-1 II- _^^^_ I III I I I I I 

-24 -22 -20 -18 -16 -14 -12 


51gft 

Figure 8. GAMA-II evolution- and density-corrected Petrosian (blue cir¬ 
cles) and Sersic (green squares) r-band LFs with best-fitting Schechter 
functions (solid lines) assuming evolution parameters for each sample as 
given in the lower half of Table[^ The dotted lines show the best-fit r-band 
Schechter functions from Table 5 of|Loveday et al.|p012|. The open dia¬ 
monds in the top panel show the ‘corrected’ LF from Fig. 7 of|Blanton et al.| 
(2005) . 


early with redshift, is able to match the GAMA LF measured in 
redshift slices. The top panel shows Petrosian r-band LFs for the 
full GAMA-II sample measured in four redshift slices as indicated, 
calculated using equation with the best-fit evolution parame¬ 
ters and radial overdensities, and taking into account the appropri¬ 
ate redshift limits. If the evolution model accurately reflects true 
evolution, and if we have successfully corrected for density vari¬ 
ations, then these LFs should be consistent where they overlap in 
luminosity. In the bottom panel we have divided each LF by the LF 
determined from the full sample (0.002 < z < 0.65; top panel 
of Fig. in order to make differences more clearly visible. We 
see that the lowest redshift LF, z < 0.1, is about 10-20 per cent 
lower than the 0.1 < 2 < 0.2 LF, indicating that the linear evo¬ 
lution model is somewhat undercorrecting at the lowest redshifts. 
The low-redshift underdensity is particularly severe at the bright 
end: the most luminous (Mr — 51g/i < —21.5 mag) galaxies are 
underdense by ~ 50 per cent relative to the higher redshift slices. 

In order to investigate these discrepancies further, we repeat 
this analysis using redshift-sliced LFs determined using Sersic 



Figure 9. Top panel: Petrosian r-band LFs measured in redshift slices as 
indicated for the full GAMA sample, applying evolutionary corrections as 
given in Table Bottom panel: The same LFs relative to the overall LF 
from the top panel of Fig.[^ 


magnitudes, with results shown in Fig.[T^ The underdensity of lu¬ 
minous, low-redshift galaxies is now much less severe; instead we 
see an increased scatter between redshifts at the bright end, with 
perhaps the 0.1 < 2 < 0.2 LF biased high relative to the oth¬ 
ers. It thus seems likely that the underdensity of luminous, low- 
redshift galaxies apparent in Fig.|^is largely due to Petrosian mag¬ 
nitudes missing a significant fraction of the flux of luminous galax¬ 
ies, which will tend to have a de Vaucouleurs-like profile. This 
problem is further exacerbated for such galaxies at low-redshift, 
which will have large angular extent, and thus also be susceptible to 
poor background subtraction: |Blanton et al.| ( |201 l| l show that galax¬ 
ies of radius ~ 100 arcsec have their magnitudes underestimated 
by around 1.5 mag in the SDSS DR7 database. Both Sersic and 
Petrosian photometry are subject to over-deblending or ‘shredding’ 
of large galaxy images. When running the Sersic fitting pipeline, 
[Kelvin et al.|p01^ aimed towards undershredding, as they were 
specifically focussed on the primary galaxies in systems with close 
neighbours. However, this does mean that the Sersic fluxes become 
susceptible to non-detection of nearby secondary sources, which 
introduces a positive flux bias in crowded fields for a small fraction 
of galaxies, see Section]^ 

In conclusion, while the redshift-sliced Petrosian LFs do 
show some systematic differences, use of GAMA-measured Sersic 
magnitudes, which capture a larger fraction of total flux for de 
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Figure 10. As Fig.|^but using Sersic magnitudes. 


Vaucouleurs-profile galaxies, and which have an improved back¬ 
ground subtraction compared with SDSS DR7, largely mitigates 
these differences, and suggests that our evolution model is a rea¬ 
sonable one. 


8 DISCUSSION 

8.1 Comprison with previous results 


While our evolution-corrected LFs agree well with previous esti¬ 
mates, our finding of positive density evolution (in the sense that 
comoving density was higher in the past) is at odds with most pre¬ 


vious work which has tended to find either mildly negative (Cool 
et al.| 201^ or insignificant ([Blanton et al.|2003[ [Moustakas et al. 


2013 density evolution. |Faber et al. j2007^ find a declining co¬ 


moving number density with redshift for their red sample, with no 
noticeable density evolution for their blue and full samples. |Zucca| 
|et al.| ^2009^ also find a declining comoving number density with 
redshift for their reddest sample; for their bluest galaxies, they find 
increasing number density with redshift. 

At least some of the discrepancy between the sign of the den¬ 
sity evolution between us and e.g. |Cool et al.| ( [2012^ might be ex¬ 
plained by the way in which the LF and evolution are fitted. |Cool| 
|et al.| ( [^12t fit the characteristic magnitude M* to each redshift 
range using the |Sandage et al.| ( |1979^ maximum-likelihood method, 
holding the faint-end slope parameter a fixed at its best-fit value for 



Figure 11. Stellar mass function for low-redshift galaxies (z < 0.06) de¬ 
termined from GAMA-II with (blue circles) and without (green squares) 
applying a weighting correction for surface-brightness and redshift incom¬ 
pleteness. Comparison data points for GAMA-I data from |Baldry et ah] 
\20\2) are shown as open diamonds. 


the lowest-redshift range. They then find the normalization <j)* us¬ 
ing the |Davis & Fluchr^ ( |1982| > minimum-variance estimator. Any 
over-estimate of luminosity evolution would lead to a correspond¬ 
ing under-estimate in density evolution, due to the assumption of 
an unchanging faint-end slope with redshift and the strong corre¬ 
lation between Schechter parameters. Although any determination 
of evolution will be affected by degeneracies between luminosity 
and density evolution, our method makes no assumption about the 
(unobserved) faint-end slope of the LF at higher redshifts. On the 
other hand, we do assume a parametric form for evolution. 

It is also plausible that the discrepancies between esti¬ 
mated evolution parameters are due to the uncertainties in 
incompleteness-correction required when analysing most galaxy 
surveys. For example, when we cap our incompleteness-correction 
weights to 5, we see a reduction in the estimated density evolution 
parameters. There are likely to be other effects leading to system¬ 
atic errors in the determination of evolution parameters, which are 
not reflected in the (statistical) error contours. 

A positive density evolution for the All galaxy sample would 
suggest a reduction in the number of galaxies with cosmic time, ei¬ 
ther through merging, or due to galaxies dropping out of the sample 
selection criteria as they passively fade. Neither scenario seems ter¬ 
ribly likely; [Robotham et al.| ( [20T4l ) see evidence for only a small 
merger rate in the GAMA sample. Perhaps a more likely explana¬ 
tion is that the apparent density evolution at low redshift is actually 
caused by a local underdensity, e.g. |Keenan et al.| (|20T3); |Whit-| 
[bourn & Shanks| ( |2014t . 

8.2 Future work 

There are several ways in which the present work can be extended. 

Having derived density-corrected Vmax values for each 
galaxy, it is then trivial to determine other distribution functions, 
such as the stellar mass and size functions, and their evolution. By 
way of a quick example, in Fig.[^we plot the stellar mass function 
for low-redshift (z < 0.06) GAMA-II galaxies, using the stellar 
mass estimates of [Taylor et al.|p011|l. In the mass regime where 
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surface-brightness completeness is high, \og{M/ Mq) -\- 2 log h > 
8 , we find excellent agreement with the earlier estimate from 
[Baldry et al.| ( |2012) using a density-defining population. The up¬ 
turn seen in the mass function helow \og{M/ Mq) -f 2 log h ~ 7 
will be sensitive to the incompleteness corrections applied; confir¬ 
mation of this feature will need to await the availability of deeper 
VLT Survey Telescope Kilo-degree Survey (VST KiDS) imaging 
in the GAMA regions. Future work will explore the evolution of 
the stellar mass function. 

The density-corrected Knax values will also be used to gener¬ 
ate the radial distributions of random points required to measure the 
clustering of flux-limited galaxy samples (Farrow et ah, in prep.) 

We plan to explore the possibility of using the [Taylor et al.| 
( |201 1^ stellar population synthesis (SPS) model fits to GAMA data 
to derive luminosity evolution parameters Qe for individual galax¬ 
ies. If the models can predict with sufficient reliability, the de¬ 
generacy in fitting for both luminosity and density evolution would 
be largely eliminated. This would also allow for the fact that galax¬ 
ies have individual evolutionary histories. 

We also plan to incorporate the environmental-dependence of 
the LF into our model. Note that the radial overdensities shown 
in Fig. are a poor estimate of the density around each galaxy 
since they are averages over the entire GAMA-II area within each 
redshift shell. [McNaught-Roberts et al.| ( |2014l l present estimates of 
the LF for galaxies in bins of density within 8/i“^Mpc spheres. 
We are currently extending density estimation to the full GAMA 
sample using a variety of density measures (Martindale et ah, in 
prep). 

This main focus of this paper has been to correct the LF and 
radial density for the effects of evolution, rather than to measure 
evolution per se. An alternative way of constraining evolution is 
to measure how the luminosity of galaxies at a fixed space density 
evolves. Via comparison with a model for the evolution of stellar 
populations (or luminosity evolution of the fundamental plane), one 
can estimate the rate of mass growth, e.g. [Brown et al.| ( |2007) . 


9 CONCLUSIONS 

We have described an implementation of the [Cole[ ( [2011) JSWML 
method used to infer the evolutionary parameters, the radial den¬ 
sity variations and the r band LF of galaxies in the GAMA-II sur¬ 
vey. For the overall population, we find that galaxies have faded 
in r-band luminosity by about 0.5 mag, and have decreased in co¬ 
moving number density by a factor of about 1.6 since z ~ 0.5, i.e. 
over the last 5 Gyr or so. When the population is divided into red 
and blue galaxies, the differences in density evolution parameter Ps 
are statistically insignificant. Luminosity evolution is significantly 
stronger for blue galaxies than for red. Evolution in the luminosity 
density evolution of blue galaxies is higher than that of red at the 
~ 5-a level. These findings are consistent with those of [Loveda ^ 
[et al.[ ( [2012^ based on GAMA-I and are as expected, since a frac¬ 
tion of galaxies that were blue in the past will have since ceased 
star formation and become red. 

While there still exists some degeneracy between the param¬ 
eters describing luminosity (Qe) and density (Pe) evolution for 
GAMA-II data, see Fig.[^ analysis of simulated data and compar¬ 
ison with a local galaxy sample from SDSS jBlanton et al.|2005) 
shows that we are able to recover the evolution-corrected LF to high 
accuracy. In detail, GAMA LFs are poorly described by Schechter 
functions, due to excess number density at both faint and bright 
luminosities, particularly for the red population. 


The density-corrected Knax values will be made available via 
the GAMA database. 
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